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Abstract 

We consider the problem of the motion of bodies in a self-gravitating 
system in two spacetime dimensions. We point out that this system can 
be mapped onto the quantum-mechanical problem of an N-body gener- 
alization of the problem of the H J molecular ion in one dimension. The 
canonical gravitational N-body formalism can be extended to include elec- 
tromagnetic charges. We derive a general algorithm for solving this prob- 
lem, and show how it reduces to known results for the 2-body and 3-body 
systems. 
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1 INTRODUCTION 



One of the oldest problems in physics is the N-body problem, which is concerned 
with describing the motion of a system of N particles interacting through spec- 
ified forces. It has found applications across a very broad range of fields, in- 
cluding astrophysics, condensed matter physics, plasma physics, nuclear physics 
and more. 

Restricted to the problem of gravitational physics it is particularly challeng- 
ing, even in the Newtonian case. An exact analytic solution is only known for 
the N = 2 case. For the general relativistic case there is no exact solution in 
three spatial dimensions even for N = 2 (although approximation techniques 
exist [I]). The main problem here is taking proper accoimt of the dissipation 
of energy via gravitational radiation, and progress here has relied on an eclectic 
blend of numerical schemes and approximation techniques. 

One of the newest outstanding problems in physics is that of quantum grav- 
ity: finding a consistent and predictive theory of relativistic gravitation that is 
fully quantum mechanical and which reduces in the appropriate (semi-)classical 
limit to general relativity. The most popular candidates to this end follow the 
approaches laid out in string theory and loop quantum gravity. However there 
is currently no fully successful resolution to the problem. 

In both cases considerable insight has been attained by reducing the number 
of spacetime dimensions. Nonrelativistic self-gravitating systems (OGS) of N 
particles in one spatial dimension have been very important in the fields of 
astrophysics and cosmology for over 30 years [2J, and 2-dimensional quantum 
gravity has been studied for over 25 years [3]. In the latter case a broad 
variety of possible theories presents itself. This is because the Einstein-Hilbert 
action is a topological invariant in two spacetime dimensions, and therefore 
has trivial field equations. It is therefore necessary to modify two-dimensional 
relativistic gravity in some way, and the most common procedure for doing so is 
to incorporate a scalar (dilaton) field. This has the effect of yielding non-trivial 
equations of motion, albeit at the price of indissoluably coupling the dilaton to 
gravity and all of the other matter fields in the system. Though this general 
procedure has antecedants from string theory, it bears a limited resemblance to 
(3-f l)-dimensional general relativity insofar as the latter has no such dilatonic 
couplings (though we note that the s-wave sector of (D -\- l)-dimensional general 
relativity is reproduced by specific dilaton gravity models [4]). Furthermore the 
nonrelativistic (c — > cxd) limits of such theories in general bear little resemblance 
to (H-l)-dimensional non-relativistic gravity [5j. 

However there is a particular way of incorporating the dilaton that success- 
fully addresses both concerns. The approach involves choosing the coupling so 
that the field equations couple the stress-energy T^^ of (non-dilatonic) matter 
to curvature in a manner analogous to that of Einstein gravity in (3+1) dimen- 
sions. Since the only measure of curvature in two spacetime dimensions is the 
Ricci scalar R, it is only the trace of the stress-energy that can act as the source 
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for curvature, suggesting the equation 

R = nT^^ (1) 

Consequently, as in (3+1) dimensions, the evolution of space-time curvature is 
governed by the matter distribution, which in turn is governed by the dynamics 
of space-time 0. Referred to as i? = T theory, it is a particular member of 
a class of dilation gravity theories on a line (or a circle, pending the choice of 
topology for an initial data set). Long regarded as a model quantum theory 
of gravity when ^ is constant [3 , what makes it particularly interesting is 
that its non-rclativistic limit is that of the Newtonian N-body system when the 
stress-energy is that of N point particles minimally coupled to gravity. Indeed, 
it can be regarded as the two-dimensional limit to general relativity [7J. 

The R = T theory therefore forms an ideal theoretical laboratory not only 
for quantum gravity, but also for studying both the OGS and its relativistic 
counterparts. In the latter case there are physical systems with dynamics closely 
approximated by the one dimensional system. For example very long-lived core- 
halo configurations are known to exist in the OGS phase space [8 . These 
are reminiscent of structures observed in globular clusters and model a dense 
massive core in near-equilibrium that is surrounded by a halo of high kinetic 
energy stars that interact only weakly with the core. Other higher-dimensional 
referents include flat parallel domain walls moving in a direction perpendicular 
to their surfaces and the dynamics of stars in a direction orthogonal to the 
plane of a highly flattened galaxy. The relativistic OGS (ROGS) has yielded 
additional insight as to how these properties are modified by general-relativistic 
effects. The two-body problem has been solved exactly in a variety of physical 
settings. The 3-body problem has been solved numerically in both the equal [9] 
and unequal |llj mass cases. Remarkably such a system shows no evidence of 
increased chaotic behaviour relative to its non-relativistic counterpart, despite 
the high degree of non-linearity in the system [TU] • 

The purpose of this paper is to explore the connection between this the- 
ory and its potential quantum- mechanical counterpart. We note in particular 
that the constraint equation of the gravity theory is identical to the stationary 
Schroedinger equation for a N-body linear atom. We show how to solve these 
constraint equations to obtain an implicit expression (called the determining 
equation) for the Hamiltonian in the N-body case, generalizing an approach 
developed in the 3-body system In the 2-body case the determining equa- 
tion can be explicitly solved since it reduces to the defining equation for the 
Lambert-W function. This yields an expression for the Hamiltonian as a func- 
tion of the separation of the particles and its conjugate momentum. In the 
3-body case the determining equation naturally suggests a generalization of the 
Lambert W- function, whose properties have been outlined elsewhere [T^]. The 
Lambert W function appears in a number of fundamental problems, such as the 
nodal structures of wave functions of atoms [13] . 

The outline of our paper is as follows. We begin by reviewing the formulation 
of the A'^-body problem in (l+l)-dimcnsional relativistic gravity. The problem 
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reduces to that of solving a single constraint equation. We then discuss how 
this equation can be mapped onto the Schroedinger equation for an A'^-body 
generalization of the molecular ion. We then provide a procedure for 
solving the constaint equation for an arbitrary number of bodies, and show how 
it reduces to previously known results for N ~ 2 and A'^ = 3. We summarize 
our work in a concluding section. 

2 Basics 

The action integral for the gravitational field coupled with N point particles in 
two spacetime dimensions is [HI HH [H] 



where R is the Ricci scalar, g^j^ and g are the metric and its determinant, Ta is 
the proper time of the a-th particle, k = SttG/c^ is the gravitational coupling, 
and ^E" is a scalar field called the dilaton. This action describes a generally 
covariant self-gravitating system (without collision terms, so that the bodies 
pass through each other), in which the scalar curvature is sourced by the point 
particles and the cosmological constant A. 

The action ^ neglects the inclusion of boundary terms, and so is only 
correct if the spacetime is compact without boundaries. When boundaries are 
present one must include both the analogue of the Gibbons-Hawking- York term 
with extrinsic curvature [16j . but also a Hamilton- Jacobi counterterm, which 
recently has been calculated from first principles [17] . Since we are interested in 
the methods required for obtaining a general solution to the field equations for 
the N-body problem, these terms are not relevant for our purposes. For a more 
complete discussion of boundary terms for the N-body problem see ref [18] . 

The field equations are obtained by varying the action with respect to the 
metric, dilation field, and particle coordinates. After some manipulation this 
gives 




R-A = kT'"" 



(3) 






where the stress-energy due to the point masses is 
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^/^=Xj"^'^ J dTa-^gfj.ag.^p^^S'^ {x - ZaiTa)) (6) 

and is conserved. Note that (|3l4p are a closed system of + 1 equations, 
which can be solved for the N degrees of freedom of the point masses and the 
single metric degree of freedom. Consistent with the conservation of T^i/, the 
left-hand side of ([6]) is divergenceless, yielding only one independent equation 
to determine the single degree of freedom of the dilaton, whose evolution is thus 
governed by the evolution of the point masses via ([5]). 
Making use of the decomposition 

V^i? = -2doiVlK) + 2di{^N^K - j~'diNo) (7) 

where the extrinsic curvature is 

K - (2iVo7)"'(2aiiVi - -/-^Nidi^ - do-f) (8) 

we can rewrite the action in the canonical form 



(9) 



where 7 = gn, No ~ (— g'"')^^, A^i — gio, and tt and 11 are conjugate momenta 
to 7 and ^' respectively. The quantities A'o and A^i are Lagrange multipliers 
that enforce the constraints R'^ = = R^, where 



RO ^ -.V77.^+2.V^.n+i^- (-^) V A^_5^ ^Pl ^ ^.J (x - (.0)) 

(10) 

R'^^t:^ -U^' + 2n' + Y, -5 {x - Za {x°)) (11) 

with the symbols ( ' ) and ( ' ) denoting da and di, respectively. 

The dynamical and gauge degrees of freedom can be identified by writing 
the generator arising from the variation of the action at the boundaries in terms 
of (^'Va/t) s-nd tt'. The coordinate conditions can be chosen to be 7 = 1 and 
n = 0; upon elimination of the constraints, the action becomes 

/ = y"rf2x|^p,i„(5(a:-Za)-H| (12) 

where 

H = dxH = dxA'if (13) 
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is the reduced Hamiltonian, with A = d'^/dx^. The field ^ = ^ {x, Za,Pa) is 
understood to be determined from the constraint equations which are now 

Avf-i^ + ^V-^ + K^V^lT^^Cx-zJ-O (14) 

a 

2Ax + J2PaH^~Za)=0 (15) 

a 

where ir — x' ■ The consistency of this canonical reduction can be demonstrated 
[14] by showing that the canonical equations of motion derived from the reduced 
Hamiltonian (|13p arc identical with the canonical field equations. 

The procedure for obtaining a solution to the system is as follows. Equations 
P4I15|) first must be solved to obtain 5* as a function of the phase-space variables 
(za^Pa), subject to the boundary conditions that the Hamiltonian (|13|) is finite 
for large |a;|. Insertion of this result into (fT5|) then yields the Hamiltonian as 
function of the physical degrees of freedom (za,Pa) of the system. Hamilton's 
equations then can be employed to solve for the evolution of the system. Using 
these results, the remaining components of the metric can be obtained by solving 
the field equations that follow from ^ [15] . 

Note that in (1 + 1) dimensions the dimensionless potentialGmr/c^ between 
a pair of particles separated by a distance r becomes infinite at spatial infin- 
ity, in contrast to (3+1) dimensions where the corresponding quantity Gm/rc^ 
vanishes at spatial infinity. Consequently the finiteness of the Hamiltonian in 
(1+1) dimensions must be ensured by an appropriate choice of boundary con- 
dition on the variables 5* and x- This condition has been shown to be [131 [TS] 

— An^x^ = for large large and guarantees both that the surface terms 
that arise in the action vanish and that the Hamiltonian is finite for the A^-body 
system. 



3 Comparison with the Schroedinger Equation 

Eq. (fT4)) is the same as the stationary Schroedinger equation for a single-electron 



atom with N field sources. Consider first the situation z„ < z„_i < 
This divides the region into n + 1 regions defined by 



< Zl. 



Region : zi < x 

Region i : z^+i < x < Zi 

Region n : x < Zn 



Setting Epi = \rpf+mf and ^E* = — 41n|0| allows us to rewrite (fT4| as 



V 



i^'^ix')^ - + KEpiS{x - Zi) 



<l){x) = 



(16) 



where the summation convention applies to the repeated Latin index. 
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Eqn. has the general solution 



X 



'-Pi\x — Zi\— eXx + eC^. 



(17) 



where X and are integration constants. The factor of e appears so that 
behaviour of x is odd under time reversal; when <: — > — t , e — > — e. In each 
region r, the function x' is spatially constant 



where we have defined 



X 



1 if i > r 
-1 ifi<r 



(18) 



Insertion of (fT5|) into eqn. (fTB]) yields 



or alternatively 



1 



A 



K (eX + -piSpt^r) - 77 + nEpi6{x ~ Zi) 



{x)^Q 



(19) 



V^f/- + - [nEpi5{x - Zi)] (j){x) ^ £ (Ep,) <j>{x) 



(20) 



(21) 



where 

A 1 

Equation (j20p is formally the same as the stationary Schroedinger equation 
for a particle of mass m and energy ^£ interacting with a linear molecule 
whose field sources have charges nEpi . This is an N-body generalization of the 
problem of the molecular ion in one dimension[Tni [201 IS]- We note that 
there exists a D-dimensional version of , which can be obtained via a scheme 



called 'dimensional scaling' [22l [23| [24| ■ The 2-body case of the one-dimensional 
limit of H J is solvable in terms of a generalized Lambert W function [15] . 

It is important to note that this analysis can be extended to include an (1+1)- 
dimensional electromagnetic field between the interacting (charged) particles |25j. 
In this case, eq. (fT6)) is generalized to [25l (46)]: 



V 



r\2 



-V 



A, 



+ KEpi5{x - Zi) 



(x) =0 



(22) 



where the potential V{x) = E^ — C, E is the electric field and Ae is the effective 
cosmological constant Ag = A — kC where C == ^a)^ where Ca is charge 

associated with the a*'* particle as defined in [551 (40)]. Thus, in the presence 
of an electromagnetic field, eqs. ([5D1) and ([5T|) are generalized to 



-V + KEpiS{x - Zi) (l){x) ^ £ {Epi) (f>{x) 



(23) 
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where 

A 1 

£ (Ep,) = —{X [Epi) + -ep.sp^^r? (24) 

o 4 4 

Note that when the gravitational couphng k is weak relative to the electromag- 
netic couplings, we recover the conventional Schrodinger equation for electro- 
magnetically charged particles. 



4 Solving the Constraint Equations 

In this section we outline the general procedure for solving the constraint equa- 
tions (fT4)) and p3|) in the A^-body case. 
Eqn. ()16|1 when x ^ Za 



V 



l\2 



V 



can be solved in each region, yielding, 
where 



(25) 



(26) 



(27) 



The solutions p6|) can be regarded as 'free particle' solutions. We procede in 
the same fashion as for a quantum mechanics problem with field sources that 
are Dirac delta functions. Outside these sources, the particles are 'free' since 
there is no potential. 

This set of equations must satisfy the matching conditions at the locations 
X = Zj. These are 



4)'^_^{Zr) ~ (i)'^{Zr) = ^KEpr(j>{Zr) 



(28) 
(29) 



From conditions 



and 



- Br-ie 2 



if. 



Br 



HEr, 



2K, 



nE„ 



2Kr 



r-l 



2Kr 



-Bre-^ 



Adding eqns pO|) and ((3T|) gives 



(32) 
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whereas 

^"■-'^ = ak;7, ak;7, ^'■^ 

(33) 

is obtained upon subtracting them. 

The boundary conditions that ensure the Hamihonian is finite in the regions 
and n is, 

- 4k2x^ + 2Ax^ = Co/„a; (34) 

where Cq and C„ are constants that have to be determined. Eqn. (I34p imme- 
diately gives 

An = So = (35) 
Using Eqn. in eqns ([5^ and with r = n then imphes 

4a„-i 4a„_i 

Defining 

^AAr = KEpr + 2(iirr + -^'^i — l) L^Br = K^/pr ^ 2(Kr — Kr-l) (38) 

LbAt — KEpr + 2{Kr - Kr-l) L b Br — I^Epr — 2{Kr + Kr-l) (39) 

CAAr = e^^^-^'--)^- eAB. = e-^(^'-+^'-)^^ (40) 
esAr = e^(^'-+^'-)- eBS. = e-^(^'-^--)- (41) 



allows us to write eqns. ([321) and (|33|) as 

^r-l = (^AAr^rGAAr + iABrSrCABr) (42) 

4Ar-l 

^r-l = { — LBArArCBAr — LsBrBrCBBr) (43) 

4Ar-l 

For r = n — 1, eqn. ([36l [37 ]) give 

= -7T7? 7? {LAA,n-lMAneAA,n~l + LAB,n-lMBneAB,n~l)Bn 

MA,n-l „ ,..s 

Bn (44) 



42i^„_iK„_2 

1 

MB,n-l 



Bn-2 — -PTF? T? { — LBA,n-lMAneBA,n-l — L BB ,n-lM Bn^^B B ,n-l)Bn 

A^Kn-lKn-2 



4^A„_lA„_2 
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Repeating this for r = n — 2 yields 
1 

— -jp ^ {LAA,n-2MA,n-ieAA,n-2 + LAB,n-2MB,n-ieAB,n-2 

Br, 





-2Kn-3 


MA,n 


-2 




-2-f^n-3 


1 








MB,n 


-2 



Bn-Z — -J^Ti^ IP {~LBA,n-2MA,n-ieBA,n-2 — L BB ,n-2MB ,n-ieBB 

43i^„_ii^„_2if„-3^" 

In this way a pattern clearly emerges, giving the general form for An-r and 
B'n—r as, 

where MA.n-r+i and MB.n-r+i are defined by the recurrence relations, 

Ma LAA,n-rMA + LAB,n-rMB .71 — 7'-\-l^AB ,n — r 

(50) 

-Ws,n-r = —LBA,n--rMA,n-r+ieBA,n-r — LbB ,n-rMB ,n-r+ieBB ,n-i{5^) 



with 

MA,n+l = MB,n+l = 1 

Evaluating the recurrence relations ([5T|) and ([50]) for the case r = n, and resort- 
ing to eqns. ([35)1 furnishes the relation 







Mb,o 



since _B„ cannot be equal to zero. 
Furthermore 



M 



B.l 







iJ = -- / dec v'' * = -- 



1 






4 


















oo 







4 1 
K 2 



K 



k\ 2 



(52) 



where the last equality holds for the case Pi\spifi\ = Pi\spi^n\ = J^Pi — 0- This 
is the centre-of-inertia coordinate system; we can always choose it without loss 
of generality. Eqn. ([5^ can always be inverted to give X in terms of the 
Hamiltonian, 

1 In^H^ A 
^ = ^«V^+2- 
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Setting X > Ohy conventiorQ and substituting back into the recurrence relations 
yields the determining equation for the Hamiltonian. 



For N ^ 2 the recurrence relations yield 



AKf + [nJpj + m\- 2Ko][kJpI + ml - 2K2] tanh -Ki\zi - Z2 



-2Ki [kJpI + ml - 2Ko\ + [t^J pi 



2K2 



(54) 



which is the determining equation for the Hamiltonian as defined via eq. ()53p . 



Note that in this particular case Kq = K2 = y t^'^X'^ — ^ in the centre-of inertia 
system. 

The = 3 case is somewhat more tedious, though the procedure is straight- 
forward. The result is 



(0 

-((' 



Ml + ifi A/3 + M: 



Ml +Ki\K. 



A/i + A/2 + Ki\[Kh + KA + Ki, K2 tanh -/C3232 



1 



((a/i + kl) (a/2 + A/3 + /q) + kj) ks tanh ( -k2Z2i 



Ml + A/2 + A/3 + Ki + KA K2K 



= 



K^ 
4 



A/3 + Ki ) iiTj ) tanh I —^2:32 I tanh 1 —-Z21 



K2 
4 



(55) 



where 



A/,; 



K\lpj + mf 



K, = -2 



\ 



X 



A 
2 



(56) 



This system for the equal mass case is discussed in detail in ref. [26|- When 
A = this equation can be solved using a generalization of the Lambert W- 
fimction [12]. 

When the Zi's are not ordered such that z„ < z„_i < . . . < zi, as can occur 
when two particles cross one another, we can treat this situation by relabelling 
the coordinates so that the ordering is in terms of decreasing magnitude. The 
algorithm just described can then be applied. 

When the z^'s are degenerate we are faced with a different situation. In 
quantum mechanics this corresponds to two charges merging into one. In the 
case of two bodies (ie. Hj) there are two solutions, gerade (symmetric) and 



*Note that solutions with X < are equivalent to solutions with X > under a parity 
transformation. 
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ungerade (antisymmetric), whereas for the H atom (the single Dirac delta func- 
tion) there is only one solution. The gerade solution maps onto the H atom 
solution continuously for both the energy and the wavefunction. However the 
ungerade solution goes into the continuum near R = I, where R is the sepa- 
ration between the two nuclei in atomic units (for more details see ref. [H]). 
It can be analytically extended below R = 1 but it corresponds to a different 
branch of the Lambert W-function. This latter energy does not map onto that 
for the one-dimensional H atom. Furthermore the ungerade symmetry means 
the wavefunction vanishes at the midpoint; as R —^ this will not reproduce 
the proper behaviour for the H atom. In one spatial dimension this problem 
can be avoided by allowing the particles to cross as discussed in refs. [TU] and 



5 Discussion 



We have shown that the constraint equations of the relativistic one-dimensional 
self-gravitating N-body problem can be solved exactly, and have described the 
procedure for carrying this out. We have given explicit results for the = 2,3 
cases. We have also shown that this system is formally equivalent to the 
stationary Schroedinger equation for a particle of mass m and energy -^S 
interacting with a linear molecule whose field sources have charges nEpi. The 
interpretation of the terms is different: in the gravitational case the special- 
relativistic energy nEpi is equivalent to the charge of the field source. 

For future work, the obvious problem to consider is quantization. Recent 
progress on this for the two-body case in the post-Newtonian limit was recently 
made, and the shift in energy levels caused by quantum-gravitational effects 
was explicitly computed [28]. In the present case, since we already know the 
Lagrangian density whose Euler-Lagrange equations are the Schroedinger wave 
equation, we can obtain a Hamiltonian density and quantize the system. This 
procedure is akin to '2nd quantization' ; it allows TRANSITIONS between 
states, with the dilaton itself acting as the agent of transition. 

Another avenue to explore is generalization of this work to higher dimensions, 
ie to the 2-1-1 case or the 3-1-1 case. Since dimensional scaling applies to 
[22l [23t [24] , it is conceivable that a scheme of dimensional scaling might apply 
to the gravity problem if dilatons are involved. Though the generalization of 
the present work has not yet been carried out for higher dimensions, we can 
nonetheless anticipate the following features. The formulation herein yields in 
effect a "decoupling" of the system of integrodifferential equations of eqs. ([3])- 
([6]) and solves for the dilation field, a result not unlike the "decoupling" of 
Maxwell's equations leading to the electromagnetic wave equation governing 
light. In view of the functional forms of the action in eq. ^ and the governing 
field equation for the dilaton of eq. (O, we can anticipate that the resulting 
field equation will have at least a Laplacian term and terms corresponding to 
an effective potential. These are the key terms for a Schrodinger equation. 
However, it is possible that extra terms, perhaps even non-linear terms, may 
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also appear. Non-linear terms could prevent a renormalizable quantum theory 
but even a non-linear Schrodinger equation - such as the equation governing 
optical solitons - is manageable. At any rate, it is quite also quite possible that 
the inclusion of an electromagnetic field, as shown here for the (1 + 1) case, may 
yield an equation linear insofar as the electromagnetic terms are concerned and 
therefore a candidate for quantization. At a minimum, a generalization to the 
(2 -|- 1) case is bound to be instructive. 

The (3-1-1) case introduces gravitons and consequently the physical situation 
changes dramatically. However, this case is by no means hopeless. It is possible 
through the use of domain walls that a (3-1-1) dilatonic formulation would also 
lead to a Schrodinger equation. If so, the (non-linear) effects of gravitational 
radiation could then be added on to a quantized Schrodinger equation, the 
situation ressembling the well-known case of the hydrogen atom which collapses 
classically but remains stable because the quantum principle prevents its energy 
from decreasing below the ground state. At this stage no conclusions can be 
firmly drawn, though it appears that a generalization to higher dimensions is 
certainly a worthwhile investigation as it promises to be feasible to at least some 
degree. 



Appendix 



Here we derive the equations of motion for the N-body system. These can 
be derived from Hamilton's equaions: 

. _ dH__ dH_dX_ _ AkX dX 
^"^ dp.] dX dpj ^J'^^x^~a9pj 

■ - - 9HdX _ AkX dX 

dzj dX dzj ^n'^X^ - - ^^■^ 

The only thing that needs to be done is to determine the partial derivatives 
of X with repect to Pa and Za- This can done by differentiating eqn. (|53p in 
conjunction with the recurrence relations. 

Before proceeding to the actual derivation, it is worthwhile cleaning up the 
notation. Differentiating (|27p with respect to Pa and Za gives 

dKr dX n^iX + lep^Sp^.r) dX ^ ._. 



dzj dzj Kr dzj 

dKr _ ( dX 1 \n'^{X + \epiSpi,r)_(dX 1 



Thus differentating the recurrence relations ([5T|) and (|50p with respect to Zj 
gives 
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dMA,n-r ^dX 

K = {L>Kr,-r + ^K^_^_i)MA,n-r+ieAA,n-r 

OZj OZj 

dX 

— '^T, {DKn-r ~ DKn-r-i)^B,n-r+ieAB,n- 

OZj 



J dMA,n-r+l 
~r J-'AA,n—r 7{ ^AA,n—r 
OZJ 



dMB,n-r+l 

+ i-'AB.n-r 7, &AB,n-r 

OZJ 

1 dX 

+ 2^^'^^^'^-'' ~ ^K„_r-i)LAA,n-rMA,n-r+ieAA,n- 

1 dX 

— {Dk„_^ + DK„_^_i)LAB,n-rMB,n-r+ieAB,n- 

A OZj 



dMB,n-r ^dX 
OZJ dzj 

dX 

+ {DK„-r + -Dj<-„_^_i)Mb „_r+iesB „_r 

OZJ 

dMA,n-r+l 

~ J-'BA,n—r 7^ ^BA,n—r 

OZj 

dMB,n-r+l 

— J-iBB,n-r 7 ^BB,n-r 

OZJ 

1 dX 

— -^-7 {Dx^-r- + DKn-r-i)^BA,n-r^A,n-r+lGBA,n- 

2 ozj 

1 dX 

+ 7:^ {Dxn-r ~ DK„-r-i)^BB,n-rMB,n-r+ieBB,n- 

2 OZj 



Repeating for pj gives 
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BMa.u- 



2 ( ^ h -T^SpjA {Dk,,_^ + DK„_^_j^)MA,n-r+ieAA,n- 



dpj \ dp J 4 

dX 1 



J dMA^7i~r+l 
+ J^AA,n—r 7\ ^AA.n—r 

opj 

dMB.n-r+1 
+ l^AB.n-r 7{ S.AB,n-r 

opj 

+ 2 ^^^^ ^ ^^^PJ,!-^ {Dk„_^ — DK^_^_i)LAA,n-rMA,n-r+ieAA,n- 

2 i^'dp ^ 4^*P-^''"^ {Dk„-,. + DK„_^_i)LAB,7i~rMB,n-r + ieAB,n- 

I ^ T TeSpj.r (^K,,-,. — ^K„-r-l)^^-'A,n-r+ieBA,n-r 



(63) 



^ / 9X 1 V 

+ 2 I + -^^^Pj.r J [Dk„_,. + DK„_^_jMB,n-r+ieBB,n-r 

dMA,n-r+l 
~ ^BA,n—r 7{ SBA,n-r 

opj 

dMB,n-r+l 
— l^BB,n-r 7{ RBB,n-r 

opj 



{DK„-r + DK„_^_i)LBA,n-rMA 



1 fdX 1 

+ 2 (^'g^ ^ ^^^PJ,!-^ {Dx^-r ~ DK,^-r-l)^BB,n-rMB,n-r+ieBB,n- 



(64) 



Combining these recurrence relations, the derivative of eqn. [53] with respect 
to zj and pj, the equations 

— = and — = 0, (65) 

OZJ Opj 

and the boundary conditions, 

MA,n+l = Mb,1 = (66) 

win give a hnear equation in and another one in , which can always be 
solved. Substituting back in eqns. (l57|) and (|58|) will give Hamilton's equations 
in terms of X. 
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